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Frogs of the species Fejervarya limnocharis and related species are 
among the most common in Southeast Asia. We studied 52 individuals from 
eight populations of the Great Sunda Islands by means of allozyme electro- 
phoresis, mtDNA sequencing and morphometry. Patterns of variation of all 
characters among populations were congruent in separating one Javanese 
population from all other populations (Java, Sumatra, Borneo) as a distinct 
species. The frequencies of four distinct mid-dorsal stripe showed no 
distinct geographic pattern. Morphometric analysis unambiguously permits 
assignment of the name Rana limnocharis Gravenhorst, 1829 to the 
widespread form that occurs in Java, Sumatra and Borneo. The two taxa 
recognised in this study occur in sympatry in Java. As the taxon known only 
from Java lacks a name, we describe and name it. 
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INTRODUCTION 


Description of existing biodiversity is one of the main goals of conservation biology. 
Especially in the tropics very little is known about the real number of species and our 
knowledge on intraspecific variability of tropical animal and plant species is even lower. The 
Amphibia of most parts of the Oriental faunal region have been poorly studied until now 
(INGER, 1999). This holds especially for the Greater Sunda Islands, although two reviews of its 
batrachofauna were recently published (MANTHEY & GROSSMANN, 1997; INGER & STUEBING, 
1997). Within the last decades, many new species have been described from these islands based 
on morphological traits (e.g., INGER, 1989; INGER & Gritis, 1983; INGER & STUEBING, 1991, 
1997). However, the few biochemical studies on Southeast Asian anurans have made apparent 
that a great amount of cryptic variation is present within species (e.g., INGER et al., 1974; 
DONNELLAN et al., 1989; Kosucu et al., 1997; Topa et al., 1998a), 


The ranid genus Fejervarya Bolkay, 1915 is distributed throughout Southeast Asia 
(Dusots & OHLER, 2000). These species show substantial morphological and colour variation 
over their distribution range. In Nepal, studies of mating calls showed the co-existence of four 
species (Dusois, 1975). In their study on genetic divergence among Southeast Asian Rana 
limnocharis, Toba et al, (1997, 1998a-b) discovered genetic divergence and sympatric cryptic 
species in this group without morphological differentiation. 


The present study focuses on the variation within frogs known as Fejervarya limnocharis 
(Gravenhorst, 1829) in the Greater Sunda Islands, a group of frogs that are known to be 
abundant in areas of intensive human activities (INGER, 1966; INGER & STUEBING, 1989, 
MANTHEY & GROSSMANN, 1997). To describe the pattern of variation in this common frog 
group we compared the within- and between-island variability of molecular and morpholog- 
ical traits. 


The nomenclatural situation in this group of frogs was presented in detail in part | of this 
series of papers (Dusois & OHLER, 2000). Generic classification follows Dusois (1992), 
slightly modified by phylogenetic data of DNA analysis (MARMAYOU et al., 2000). In order to 
formalise the systematic position of the different forms of Fejervarya found in our genetic 
analysis, we linked molecular data to morphometric data. This allows us to cross-reference old 
names based on museum specimens with results of modern molecular population studies. 


MATERIAL AND METHODS 


Frogs were collected from six natural populations or bought from local fish markets of 
Java (J), Sumatra (S) and Kalimantan (K), the latter being the Indonesian part of Borneo (fig. 
L): JI (4 specimens, FMNH 256721-256724), swampy meadow along the road between Bogor 
and Parung: J2 (11 specimens, FMNH 256725-256733, MNHN 1997.4916, ZMFK 68867), 
paddy field at Chianjur; KI (2 specimens, FMNH 256734-256735), fish market at Pontianak: 
K3 (12 specimens, FMNH 256736-256747), paddy field at the village Desa Lape; Si (13 
specimens, FMNH 256749-256761 ), fish market at Medan; S2 (2 specimens, FMNH 256762- 
256763), paddy field and small brook at village Sidikalang: $5 {5 specimens, FMNH 256764- 


Source : MNHN, Paris 


VEITH, KOSUCH, ORLER & DuBois El 


Borneo 


Java 


Fig. 1. - Sample localities of Fejervarya on the Great Sunda Islands, Indonesia (* refers to the sample 
locality 3/4 of Topa et al., 19984). 


256768), frogs collected by inhabitants at a small village about 10 km SE Tapaktuan; $6 (3 
specimens, FMNH 256769-256771), frogs collected by inhabitants of the village Desa Seleu- 
kat near Tapaktuan. In the molecular analyses we used the ranids Fejervarya cancrivora 
(Gravenhorst, 1829) (Kalimantan, Borneo), Rana temporaria Linnaeus, 1758 (Germany) and 
Rana catesbeiana Shaw, 1802 (from a frog farm at Java) as outgroups. 


A specimen from the collection of Heinrich Kuhl from Java (RMNH 4287), designated 
as neotype of Rana limnocharis Gravenhorst, 1829 by Dunois & OHLER (2000), was included 
in a discriminant analysis for nomenclatural decisions. 

Sex of all specimens was determined either by presence of secondary sexual characters or 
by observation of gonads through small lateral or ventral incision on the frogs (most of the 
frogs were ventrally opened because of tissue dissection in the field), 


ALLOZYME ELECTROPHORESIS 
Pieces of fresh muscle and liver were dissected and stored in liquid nitrogen. Frogs were 


preserved for determination and morphometry in 70% alcohol after pre-fixation in 10 Yo 
formaldehyde. 
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Allozyme electrophoresis was carried out on muscle and liver homogenates on cellulose 
acetate gels (RICHARDSON et al., 1986). Twelve enzyme systems providing information on 15 
presumptive gene loci were stained according to HEBERT & BEATON (1986). Gels were run in 
three different buffer systems: tris-glycine, pH 8.5 [acon- and -2 (E.C, 4.1.1.3), fum (E.C. 
4.2.1.2), aat-l and -2 (E.C. 2.6.1.1), mdh (E.C. 1.1.1.37), PePay-Leu (E.C. 3.4.11), pgm (E.C. 
5.4.2.2) and tpi (E.C. 5.3.1.1)], tris-citrate, pH 7.0 [idh-1 and -2 (E.C. 1.1.1.42), me (E.C. 
1.1.1.40) and mpi (E.C. 5.3.1.8)] and phosphate, pH 8.0 [/dh (E.C. 1.1.1.27) and pk (E.C. 
2.7.1.40)]. Average expected heterozygosity (H,) and mean number of polymorphic loci (P) 
were calculated as population characteristics using the computer program G-STAT (Siais- 
MUND, 1993). NEt's (1972) standard genetic distance estimates and UPGMA cluster analysis 
were calculated using the program NTSYS (RoutF, 1990). 


DNA-SEQUENCING 


DNA was extracted from frozen muscle tissue using the standard phenol-chloroform 
protocol of SAMBROOK et al. (1989), Double-stranded PCR amplification was performed in 
50 xl reactions containing | unit of Taq polymerase, 5 pl of 10 x reaction buffer (Boehringer), 
250 moi each of dGTP, dATP, dTTP and dCTP, 20 pmol each of light- and heavy-strand 
primers and 1 y] of mtDNA extract. The primers 16Sa-L and 16Sb-H (KocHER et al., 1989) 
amplified a double-stranded DNA segment of 560 bp of the mitochondrial 168 ribosomal 
RNA gene. Thermal cycling was carried out in a programmable heating block (Perkin Elmer 
Gene Amp PCR System 9600) for 35 cycles (initial denaturation step of 90 s at 94°C, 
denaturation for 45 sat 94°C, primer annealing for 45 sat 55°C and extension for 90 s at 72°C). 

PCR products were purified with a PCR purification kit (QuiaQuickspin) and directly 
sequenced with a sequenase kit (Amersham) using only primer 16Sa-L. The products of the 
sequencing reaction were resolved by the automatic sequencer 377 of Applied Biosystems. 


After sequence alignment with CLUSTAL-W (HIGGINS & SHARP, 1993), the alignment 
was adjusted manually. The aligned 390 bp sequence corresponds to nucleotides 4039-4429 in 
the Xenopus laevis mitochondrial genome (Rog et al., 1985). Phylogenetic relationships were 
determined using maximum parsimony (PAUP, version 3.1.1; SworrorD, 1993) and maxi- 
mum likelihood approaches (PHYLIP, version 3.5c; FELSENSTEIN, 1993). For maximum 
parsimony analysis, 1000 bootstrap replicates were run to test for confidence in the topology 
of the phylogenetic trees (FELSENSTEIN, 1985). 


Sequences have the GenBank accession numbers AF3468 10-AF346811, and the EMBL 
accession numbers AJ292014-AJ292023. 


MORPHOMETRY 


Twenty-seven measurements of 39 adult specimens were taken with a slide caliper to 
0.1 mm precision or for measurements smaller than $ mm with an ocular micrometer to the 
nearest 0.01 mm (app. 1). To control for isometry, all measurements were divided by snout- 
vent length (SVL), expressed as per thousands of SVL, or transformed into the Neperian 
logarithm. Male and female specimens were grouped after their homogeneity was established 
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(Outer, 1996). For morphometrical analyses we used SPSS statistical programs for personal 
computers (Norusis, 1992) 


Morphological analysis should test if genetically close populations from different islands 
are also morphologically homogenous and, conversely, if the two genetically differentiated 
species can be distinguished by their morphology. Non-parametric Kruskal Wallis test was 
performed to measure morphological variation in Fejervarya specimens of different islands. 
Differentiation between genetically distinct samples was tested using non-parametric Mann- 
Whitney U test. 


Discriminant analysis 


Groups obtained by the molecular analyses were subjected to a disermmmant analysis 
(Norusis, 1992) The measurement data for the female neotype of Rana limnocharis were 
incorporated into the analysis after the production of the original discriminant function 
model based on all other data, following examples of nomenclatural clarification using this 
approach (HEYER, 1994; OHLER, 1999; OHLER & Dusots, 1999). 


Colour patterns 


Frogs of the genus Feyervarya are known to frequently have mid-dorsal stripes, which are 
an mteresting character for evolutionary biology studies (sce eg MILSTEAD et al. 1974; 
Dubois, 1980). Heritability studies by MOoRIWAKI (1953) and Monanty & Dutta (1999) 
suggested that presence or absence of a dorsal line in Feyervarya is coded by two alleles of 
a single locus. The allele for the striped phenotype appears to be dominant over that for 
the unstnped one. However, Moriwaki (1953) and others (e.g, SHIBATA, 1988) only 
discriminated between two phenotypes (striped and unstriped). Dusots (1977) and MOHANTY 
& DUTTA (1999) distinguished between three phenotypes. a pattern showing no line at all, a 
pattern showing a fine line and a third pattern showing a wide stripe The patterns are not 
distributed equally m the four different species of Fejervarya from Nepal which can be 
distinguished by their mating calls. whereas F teratensts (Dubois, 1984) and F pierrer 
(Dubois, 1975) show all three different patterns, F sifadrensrs (Annandale, 1919) does not 
show the wide stripe phenotype, and all F nepalensis (Dubois, 1975) specimens have mid 
dorsal stripes (Dusois, 1977) 

Four mid-dorsal patterns can be recognised m the Feyervurya specimens that we studied 
(fig 2} (1) no mid-dorsal line or stripe. (2) a fine clear mid-dorsal line alone present, (3) a wide 
mid-dorsal stripe alone present, (4) both a fine hne and a wide stripe present, superimposed. 
Presence and absence of both kinds of hnes as well as frequency of combination of the two 
lines and their distribution among populations were noted, and their pattern of variation was 
studied using Kruskal-Wallis test. 
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Fig 2. (a) Holotype of Fesertartecishandart MNHN 1997 4916, from Chianjur, Java, phenotype no hne/no stripe, 168 haplotype J2-C, (b) Feyervarya 
tumiocnaris from Desa Lape. Borneo, phenotype hne/stripe, 16S haplotype K3 A. (c) Fezerrarva limnocharis from Tapaktuan 2, Sumatra, phenotype 
Ime/no stripe, 165 haplotype S5-A. Drawings by Kathe Rehbinder, Mainz. 
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RESULTS 


ALLOZYME VARIATION 


Of the 15 allozyme loci studied, five were monomorphic (fum, rdit-1, 1dh-2, me, pk) within 
and between the Fejervarya samples. Acon-1, acon-2 and Idh were monomorphic within 
populations, but exhibited a fixed allelic difference in population J2 relative to all other 
populations (app. 2) Between populations J1 and J2 further allehc differences appeared to be 
fixed at the loc: aat-2, mp: and pep. Thus, these two geographically closely related populations 
did not share alleles at six out of 15 loci. 


Intra-populational genetic variability (47,) was lowest in the Sumatran samples with the 
exception of sample SI from the Medan fish market. This was the only sample that showed 
significant genotype deviations from Hardy-Weinberg expectation at two loci (mpr and pep; P 
< 0.05) due toa lack of heterozygous specimens, indicating a mix-up of different populations. 
Sample $5 was completely monomorphic Three to six loci were polymorphic in the Javanese 
and Borneo samples, even though sample sizes were in some cases lower than for S5. 


Net's (1972) standard genetic distance between populations ranged from 0.004 between 
Sl and $2 and 0.337 between J2 and SS (tab. 1) Population J2 was the genetically most 
distinct. It was separated even from the geographically close population JE by a relatively high 
genetic distance estimate of 0.316. The Bornco and Sumatra populations formed distinct 
clusters in the UPGMA phenogram (fig 3). 


MIDNA SEQUENCE DIVERGENCE 


We sequenced 390 bp of the mitochondrial 16S rRNA gene from 52 mdividuals of 
Fejervarva and from the three outgroup species. One hundred twenty seven bp pos.tions were 
variable. Samples of Feyervars a comprised nine different haploty pes (tab 2) with 60 sites being 
variable among them, Fifty-two of these variable sites were parsimonious informative, 49 of 
which supported the monophyly of either of the two main haplotype lineages. 


The average substitution rate between the haplotypes of population J2 and the remaining 
populations from Java, Sumatra and Borneo was 0135 This 18 only about 0 02 less than 
between either of the two Fejervaria forms and F cancrnora, thus suggesting the differenti- 
ation of J2 from the other Fejorvarra fimnochia ts populations as distinct species. 

The maximum parsimony and maximum hhelhood trees produced exactly the same 
topology with two rather homogeneous but well differentiated groups of Fejerrar a mito- 
chondrial 16S haplotypes (tab 3) Bootstrap support in the maximum parsimony tree was 
100 for both clades (fig 4) which we subsequently name Ji and J2 No further significant 
topology was found within the two clades. 
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Table 1 - NEP's (1972) standard genetic distance estimates between 9 samples of Fejervarya 
from Java (J), Kalimantan (K) and Sumatra (S) See text (Material and methods) for 
details on localities Sample Kc belongs to the species Fejervarya cancrivora 


OI CO fa fe pa fe ES EI 


0188 | 0282 


0169| 0298 


0.102] 0319 


0 182 | 0324 


Table 2 - Variable bp positions of a partial sequence of the mitochondrial 16S-rRNA gene of nine 
Fejervarya and three outgroup taxa haplotypes. 


Fe arvarya 
Fejervarya 
Fejervarya 
Fejervarya 
Pejervarya 
Fejervarya 
Fejervarya 
Fejervarya 
Fejervarya 


F. cancrivora Ke T..C.TRA. 


J2B 
J2-A 
J2-c 
Ji-A 
31-8 
K1-A 
K3-A 
S5-A 
S6-A 


Rana temporaria 


90000000000000000013121111111221111122222222222222222222222222222222 
112445556667788898004446888889999999000000001111111112222222222333 33 
580074673586725670141574145 781346 789012 34674012 146 7890123456 78912 156 


CGCATAGGACGCCACCTCTOGTACCATAGACCAATCAACAC ACTATOTAATROTTICTOC COCGCTTE 
+ + + T.GT....TAA.....GACT.G.CT,.TCT..C.ATACC.. AAGIT TOA 

+ T GT ...TAA.....GACT G.CT .TCT..C AAACC AAGTT + TCA 

++ APOT TAN GACT.G.CT..TCT.,C.AAACC.. AAGTT TCA 
Bec? TAN GACT.G.CT..TCT..C.AAACC.. .AAGTT - TCA 

eb TGT TAA GACT G.CT .TCT .C.AAACC AAGTI "TCA 

. +T.GT,..,TAA.... GACT G CT TCT CARAC AAGTT -- TCA 
«TG... TC, T.C...CGA..TT.G CC TC CCCAA C AGG TOGO 


-A.C.. ACTC, 1G, TAT. TCA.GICCC.CT.GTG.CC-T...CC.TA.CACAAGAGAT A LC 


Rana catesbeiana — .A.CC.A.TIC..G.TATC.T...TTCTACT..TG CC-T- TA. . A. . TCACAAGAAR TAGC 
4224222242222222222221311VI4V131MI34341II344413IsALoIADALE 


Fejervarya 
Fejervarya 
Fejervarya 
Fejervarya 
Fejervarya 
Fejervarya 
Fejervarya 
Fejervarya 
Fejervarya 


F, cancrivora Ke 


J2 B 
J2-A 
ac 
JI-A 
J-B 
Ki-A 
K3 A 
S5-A 
56-A 


Rana temporaria 
Rana catesbeiana AG.TITAA.TA.AAT CA.C..T.C.ATAC..TACTC A T AA C AT 


35566777788888899999900001111112233333334444444445555566678 
83869013912456823567807892347895901267890123456781478901800 


GTAACATT-CGCGTCA-TA-ACCAAGCCCGCGCGACTTGACGTTTTTTTTAATATTCAC 
CA Y CCCT MA ATT A ca a 
CNET C.C.CT.TTA..ATT. A. ČA To 
CENITI A cA A 
AT «CC CT.TTA ATITA GA y 
A, CCCTTIA ATT A sa ot 
CT.A..T......C.C.CT.TTA..ATT..A OR G 
AGE ESA A Te an Dia AA TTC 
AG.T T.ACTA.AAT-CA C....T.A.A...TACT O al 


Source . MNHN, Paris 


VEITH, KosucH, OHLER £ DUBOIS 13 


Table 3 — Average substitution rate between samples of the two Fejervarya lineages JI and 
J2 and three outgroup taxa, transitions and transversions were weighted equally, gaps 
treated as fifth base 


= [iran [iran 


JI lineage (J/K/S) | 0.006 +0004 
J2 lineage 0.135 + 0.004 | 0.006 + 0.003 


F. cancrivora | 0 160 +0 001 | 0155 +0 001 


R. temporaria 0219+0002 | 0229+0001 


R catesbeiana 0 209 + 0 002 | 0212 +0002 


MORPHOMFTRIC VARIATION 


For all measurements and ratios, the mean, standard deviation and minimum and 
maximum values are given for adult specimens of Feyersarya separating samples by sex and 
geographical origin (app. 3) Comparisons of adult individuals using the Kruskal-Walhs test 
(tab. 4) show significant differences in the frogs from the three islands in the ratos which all 
concern head shape (head length, distance between anterior border of eyes, distance between 
posterior border of eyes, internartal distance) Fejervarya from Java have significantly shorter 
heads and more pointed snouts than the frogs from Sumatra and Borneo. Frogs from Sumatra 
are significantly distinct in having a greater distance between posterior border of eyes (app. 3, 
tab. 4). 


Dorsal pattern 


In the populations we studied, 32 1 % of the specimens have no mid-dorsal stripe, a fine 
mud-dorsal line ts present m 57.6 ° o of all specimens examined, a wide stripe in 30 5% of the 
specimens, and the combined fine-wide stnped phenotype was observed in 153° of the 
specimens (tab. 5). 

More males lack a mid-dorsal Ime and more females have a combmation of both 
phenotypes (chi-square test, y? — 23 0181, df = 3, P < 0 001). Neither mates nor females differ 
in a significant manner from the overall distribution of dorsal pattern (chi-square test: males, 
4 5.7635, df = 3, P > 0 05; females, y'=5 1466, df = 3, P > 0 05). 

The different populations can be put m three groups a series of populations showing no 
wide stripe phenoty pes (group | in table 5}, a population from Sumatra showing no narrow 
line pattern (group 3 m table 5). and a group of populations showing all four patterns (group 
2m table 5) Statistical test shows heterogeneity of phenotype distribution for the wide stripe 
pattern between populations studied (Kruskal-Wallis test. y' = 22 38. df = 7. P < 001) 
Distribution between populations for narrow line pattern 1s not statistically heterogeneous 
(Kruskal-Wallis test: y” = 1205. df = 7. P > 005} The population from Chiangur ts m the 
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Nei’s (1972) standard genetic distance 
032 0.24 0.16 008 000 
H + + 


J2 


Fig 3 - UPGMA phenogram of Ni t's (1972) standard genetic distances between the Fejervarya samples 
(allozyme data), the tree 15 rooted by the outgroup Fejervurya cancrnora, the cophenetic value 15 0 99 


same dorsal pattern group as the population from Desa Lape (Borneo) and the population 
from Medan (Sumatra) The two Javanese populations are not in the same dorsal pattern 
group. 

In conclusion, populations of Fejervarya from the Sunda Islands show statistically 
significant variation of dorsal colour patterns, but this variation is not congruent with 
allozyme or mtDNA variation as shown in this study. 


To which population does the name Rana limnocharis Gravenhorst, 1529 appl? 


The name-bearing type provides the objective standard of reference by which the 
application of the name it bears 1s determined {Anony mous, 1999: Article 61.a) Generally, 
taxa have been described using morphological methods, and new techniques are rarely 
applicable to dead museum specimens. Systematics 1s, however, a largely historical science, as 
new results are added to previous knowledge One scope of modern science therefore is to 
define methods to find quantifiable links between historical museum specimens (especially 
type-specimens) and new material studied by modern methods. 


Muluvariate statisucs have been used to refer single type specimens to biologically 
defined populations (Htvir, 1994, OHLER, 1996. 1999, OHLER & Dusors, 1999) This also 
applies to the case of the two Fejes varva lineages here defined by allozyme and DNA data. 


The historical specimen of Hemneh Kuhl, who first collected the species Rara fusmoc ha- 
ris, was a Judicious choice for a neotype (see Dusors & OHLER, 2000) In order to find out to 
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Table 4 - Comparison by Kruskal Wallis test of snout-vent length (SVL) and of ratios of measurements m adult 
specimens of Fejervanya from different ongns. For cach sample, munumum and maximum values, mean 
and standard deviation are given. df, degree of frecdom, n sample sue, P, probability, *. significance 


level P< 0 05 


FL/SVL 


FN/SVL 


FLL SVL 


FOL) SVL. 


FTL/SVL 


HAL / SVL. 


ML/SVL 


HW/SVL 


IBE/SVL 


IFE /SVL 


IMP /SVL 


IN/SVI. 


TL SVL 


MRB SVL. 


MEE /SVL 


MN/SVL 


SVL 


TL SNL 


TL/SYI, 


TYD SVL 


Pre su 


WHE SA, 


WI SVL 


WILISYL 


WIE SVL 


96-117 
1082 106 
86-91 
884255 
196-221 
2094 126 
496-534 
5234292 
303-324 
H62116 
176-221 
1994225 
329-378 
356+241 
342-300 
30-934 
210 246 
212+188 
129-145 
1384860 
44-56 
502626 
66 ~ 69 
CHE 
96-126 
1082 154 
140-144 
145+81 
221-255 
2974175 
106-320 
114107 
444552 
4844591 
114-130 
1232797 
489-532 
042239 
63-66 
652186 
45-40 
372245 
73-102 
211165 
75-87 
82:39 
55-79 
644127 


86 102 
96:112 


97-129 
1124859 
30 101 
914545 
191-227 
2151912 
490-583 
132204 
292-338 
319+ 144 
173-219 
1972132 
369-448 
3882178 
342-365 
3472 103 


208-249 
254118 


135 ~ 162 
1472887 
45-57 
s1232) 
68-79 
734313 
103-123 
164495 
137-171 
498110 
226-294 
2492160 
310-393 
3302 192 
326-590 
48638 30 
192-130 
118467 
475-537 
514 171 
58-69 
ann 
32-44 
394332 
63-92 
772826 
65-88 
712739 
40-75 
60 +897 
74-100 
87 £869 


109-126 
152400 
81-95 
894470 
207-230 
2212680 
465-560 
$30 2241 
27-44 
192163 
180-223 
2032120 
63-416 
3914145 
319-367 
3484593 
27-271 
2394148 
141 17 
156 101 
46-57 
512349 
66 88 
154469 
108-125 
1182538 
135-193 
1584173 
229-279 
2522147 
204 373 
154 162 
162-542 
4011606 
105-145 
123 74 
465-514 
508 4169 
59-73 
66+40< 
27-49 
40474 
52-97 
782122 
61-97 
794878 
46 79 
60.942 
77 104 
CORTE 


Kruskall-Walls test 
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F. cancrivora 


Lo _____ R. temporana 
R catesbeiana 


Fig 4 Maximum parsimony tree of 390 bp of the mitochondrial 16S rRNA gene of nine Fejervarva 
haplotypes. Tree length is 202 steps (one shortest tree only), bootstrap values > 50% for 1000 rephcates 
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Fig $ Stacked histogram of diserimmant function for specamens of population J2 and other popula- 
toms of Fejervaria trom Sunda blands The neotype RMNH 4287 of Rana finimiharis Gravenhorst. 


1829 was included into this analysis without a priori group membership 
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Table 5 Coloration pattern of dorsum in Fejervarya from Sunda Islands Frequency distribution 
of phenotypes of mid-dorsal line and alleles supposed present m populations studied n, 
allele wild narrow line, N, allele narrow line, w, allele wild wide Ime, W, allele wide hne 
Proposed phenotypical groups 1, line absent or narrow line present (wide line never 
observed), 2, hne absent or fine and or wide line present (all 4 phenotypes observed), 3, line 
absent or wide line present (fine hne never observed) 


Alleles 
Fmeand | inferred to | Proposed 
wide tne | be present | group 

in genotype 


Sample (sample size) | Linc absent 


J1 Bogor (4) 


82 Sidikalang (2) 


S5 Tapaktuan (5) 


S6 Desa Seleutkat (3) 2 


S1 Medan (13) 1 


32 Chianjur (11) 5 


K3 Desa Lape (13) 


Ki Pontianak (2) 


which genetically defined lineage (population J2 or the other populations from Sunda Islands) 
this name applies, a discrimmant analysis was performed and subsequently applied to the 
neotype of Rana limnocharis. The analysis clearly classes the neotype with lineage Jt and not 
with the specimens from population J2 (tab. 6, fig. 5). This indicates that the name Rana 
limnocharis Gravenhorst, 1829 should apply to the widely distributed taxon occurring m Java, 
Sumatra and Borneo, and that a new name should be created for the species from population 
J2. As no other name ts available for this taxon (see DU Bois & OHLER, 2000), we describe it as 
a new species. The morphological comparison of adult males of the new species and 
Fejervarya lumnocharis is summarised in table 7. 


Fejervarya iskandari sp. nov. 


Diugnosis Medium sized Fejervarva with a relatively wide head and mterorbital distance. 
small tympanum, short eye length, short forcarms, short inner toes and small mner metatarsal 
tubercles. 
Description of holotipe MNHN 1997 4916, adult male (fig 2a), from paddy field at 
Chianjur (Java) (6*12'S, 107*08'E). 

(A) Size and general aspect (1) Specimen of rather small size (SVL 40.4 mm), body 
moderately stout. 
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Table 6 — Results of canonical discriminant analysis between genetically determinated 
specimens of Fejervarya 


A Statistical mgnifieance 


im SS mom ko MON 
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Neotype of Rana lummochars (ungrouped) 
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Table 7 Comparison of adult males of Feyervarya iskandari sp. nov and Fejervarya 
lmmocharıs (Gravenhorst, 1829) by Mann-Whitney U test a, sample size, P, 
probability, U, Mann-Whitney U, ns, significance level P > 0 05, **, significance level 
P<001 


Ten 
P= 0 0045 +* 
U= 
P= 0 5698 ns 
2212436 U=0 
218-230 P= 00045 
52422716 | 53641262 U=14 
489-557 522-557 | P=05698 ns 
31311344 3245837 U= 16 
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224874. | 20341574 UN 
200 224 179-223 | P=02912m 
23. 
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612287 651534 t= 
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(B) Head. (2) Head of medium size, longer (HL 15.4 mm) than wide (HW 14.4 mm: 
MN 13 1 mm; MFE 9.5 mm; MBE 5 9 mm), convex. (3) Snout pointed, protruding, its length 
(SL 7.90 mm) much longer than horizontal diameter of eye (EL 4.34 mm) (4) Canthus 
rostralis rounded, loreal region very concave, vertical. (5) Interorbital space flat, smaller (IUE 
2 27 mm) than upper eyelid (UEW 3 57 mm) and internarial distance (IN 2.92 mm); distance 
between front of eyes (IFE 5.7 mm) about half of distance between back of eyes (IBE 
9 5 mm) (6) Nostrils oval with flap of skin laterally, closer to up of snout (NS 3 55 mm) than 
to eye (EN 4 02 mm). (7) Pupil rounded. (8) Tympanum (TYD 2.72 mm) distinct, rounded, 
more than half of eye diameter: tympanum-eye distance (TYE 1.30 mm) about half its 
diameter (9) Pineal ocellus present, between anterior border of eyes. (10) Vomerine ridge 
present, bearing few small teeth, between choanae, with an angle of 45° to body axis. less close 
to choanae than from each other, longer than distance between them. (11) Tongue large, 
rounded, emargimate (12) Supratympantc fold distinct, from eye to shoulder. (13) Parotoid 
glands absent, (14) Cephalic ridges absent. (15) Co-ossificd skin absent, 


(C) Forehmbs. - (16) Arm short, rather thin (FLL 8.7 mm), as long as hand (HAL 
8.7 mm), not enlarged (17) Fingers I, IN and IV long, finger TI short, all thin (TFL 5 32 mm) 
(18) Relative length of fingers, shortest to longest: 11 < IV < I < IIL. (19) Tips of fingers 
pointed. (20) Finger II with dermal fringe; webbing absent. (21) Subarticular tubercles 
prominent, rounded, single, all present (22) Prepollex oval, prominent; two oval, distinct 
palmar tubercles; supernumerary tubercles absent. 

(D) Hindlimbs. - (23) Shank three times longer (TL 20 6 mm) than wide (TW 6 0 mm), 
longer than thigh (FL 18 8 mm), but shorter than distance from base of internal metatarsal 
tubercle to tip of toe IV (FOL 22 5 mm) (24) Toes long, thin: toe IV long (FTL 13.6 mm), more 
than one third of distance from base of tarsus to tip of toe IV (TFOL 29.6 mm). (25) Relative 
length of toes, shortest to longest’ I< I < V < TII < IV (26) Tips of toes pointed (27) Webbing 
moderate 11 2111-21111] -22/31V 21/3- 1 V(WTF 4.28 mm: WFF 3.95 mm; WI 3 89 mm; 
WII 3 24mm, MTTF 9.7 mm; MTFF 10 3 mm; TFTF 9 5mm; FFTF 10 0 mm) (28) Dermal 
fringe along toe V present, from tip of toe to base of metatarsus, well developed, (79) 
Subarticular tubercles prominent, oval, simple, all present (30) Inner metatarsal tubercle 
long, very prominent, its length (IMT 1.94 mm) less than 2.5 times ın length of toe L (TL 
4.60 mm), (31) Inner tarsal ridge well developed, along distal half of tarsus. (32) Outer 
metatarsal tubercle prominent, supernumerary tubercles absent, tarsal tubercle absent, 

(E) Skin. - (33) Dorsal and lateral parts of head and body. snout smooth, between the 
eyes and side of head with few, flat glandular warts; back with glandular folds and glandular 
warts between them. flanks with glandular warts. (34) Latero-dorsal folds absent (35) Dorsal 
parts of limbs forelimbs smooth; shank, thigh and tarsus with glandular warts. (36) Ventral 
parts of head, body and limbs: throat, chest and belly smooth, thigh with glandular warts. (37) 
No macroglands. 

(F) Coloration in alcohol (38) Dorsal and lateral parts of head and body brown with 
blackish paired spots. shoulder spots mdistinet, four brown spots on each side of upper lip 
(39) Dorsal parts of limbs forelimb, thigh, shank and foot light brown with dark bands. 
posterior part of thigh with brown and whitish marbling (40) Ventral parts of head, body and 
limbs chest, belly and thigh cream white: throat cream white with grey W-shaped pattern. 
margin of throat white with large brown spots. 
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(G) Male sexual characters. ~ (41) Nuptial spines in one single patch on prepollex and 
finger I; numerous, indistinct, cream coloured spines. (42) Vocal sacs as greyish, folded skin on 
the two sides of the throat; slit-hke openings in rather anterior part of mouth floor. (43) Fine 
horny spinules on the anterior border of the throat. 


(H) Part of the 16S rRNA sequence that corresponds to the nucleotides 3994-4554 of 
Xenopus laevis (Roe et al, 1985) (44) (EMBL AJ292018) TCTTGTTTTTTCATAA 
GAGGTCCAGCCTGCCCAGTGACACAATTAAAGGCCGCGGTACCCTGACCGTG 
CGAAGGTAGCATAATAACTTGTTCTTTAAATGGGGACTAGCATGAACGGCAC 
CACGAAGGCCTCACTGTCTCCTTTTTCCAATCAGTGAAACTGATCTCCCCGTG 
AAGAAGCGGGGATGATAATATAAGACGAGAAGACCCCATGGAGCTITAAACC 
CAATAAGCAACCCTAATCAACACAACTCATCTAAATTCTTTCTCCCCCTGCTTT 
TIGGTTTTAGGTTGGGGTGACCACGGAG IAAAACATATCCTCCACGACGTAC 
GGATTAACCCTTATCCAAGAGCCACCGCTCTAAGAATCGACAAATIGACGITT 
TTIGATCCAATATATTGATCAACGAACCAAGTTACCCTGGGGATAACAGCGCA 
ATCCATTTTTAGAGCCCCTATCGACGAATGGGTTTACGACCTCGATGTTGGAT 
CAGGGTACCCAAGTGGTGCAGCCGCTACTAATGGTTTGTTIGTTCAACAATT 
AAAACCCTACGTGATCTGAGTTC. 


Paratopotypes. FMNH 256725-256733, 3 adult males, 1 adult female, 5 juvenile males; 
ZFMK 68867, adult male. Same collection data. 


Etymology, - The new species 1s dedicated to Djoko Iskandar, herpetologist from Indonesia 


Discussion 


Fejervarva iskandart sp. nov. and Fejemarya lunnocharts (Gravenhorst, 1829) from Java 
clearly represent two genctically distinct lineages. Their Ne: D values calculated from allo- 
zyme data (average 0.314) falls within the interspectfic range known from other ranids (c.g , 
Case, 1978, NISHIOKA & SUMIDA, 1990; Mens et al., 1992; Arano et al , 1993; Brier. et al, 
1996, ViitH, 1996) However, the observed morphological variation between the genetically 
well differentiated lineages F phandarr and F lunnocharts ts shghtly more pronounced than 
among members of the other populations, as estimated by the number of significantly distinct 
measurements. Therefore these two forms are clearly a patr of sibling species that are 
morphologically scarcely distinct from one another, but show substantial genetic differentia- 
tion In addition, a comparison of our allozyme data with those of the samples 3 and 4 of 
TODA et al. (1998u) reveals that their samples fit almost perfectly with our analysis in nine out 
of eleven loci which were analysed in both studies the results are identical, one locus shows 
similar allele frequencies and only in one locus the results are different Discrepancies m the 
results may have been caused by using different protocols m the two laboratories (sec Vita, 
1994 for further examples) Therefore the new species includes specimens from the Javanese 
populations of Chianjur and Malingping (sample 4 in Toba et al . 19984). Its genetically well 
defined, but morphologically very slar to £ Anmocharis Distribution ranges of F ihan- 
dartand F lanodharn overlap at least between Malingping and Chianjur over a distance of 
va, (30 km Syntopy was shown for a paddy field near Malingping, Java (Tova et al , 199840), 
confirming species status of F iskanduri. 
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Discriminant analysis of morphometric data resolves the nomenclatural relevance of the 
observed genetic divergence The name Runa limnocharis Gravenhorst, 1829 applies to a 
taxon widely distributed in the Sunda Islands. As there is no indication of genetic or 
morphological differentiation of the frogs from Borneo island, the name Rana was! Annan- 
dale, 1917 (based on a holotype from Kuching, Sarawak, Malaysia, on this island) 1s here 
considered a junior subjective synonym of Rana Imnocharis (see DuBois & OHLER, 2000), but 
remains available for further systematic decisions, should Bornean Feyervarya prove hetero- 
geneous 


RÉSUMÉ 


Les grenowiles de l'espèce Fejervarya limnocharis et d'espèces voisines de celle-c1 sont 
parmi les plus communes en Asie du Sud-Est. Nous avons étudié 52 spécimens de huit 
populations des Iles de la Sonde par les méthodes de l'électrophorèse d'allozymes, du 
séquençage d'ADN mitochondrial et de la morphometrie Les patrons de variation de tous les 
caractères entre les populations se montrent congruents pour séparer une population de Java 
de toutes les autres populations (Java, Sumatra, Bornéo) et considérer qu'elle appartient à une 
espèce distincte Les fréquences des quatre phénotypes observes concernant la présence et la 
largeur d'une ligne médio-dorsale dans les populations exanunécs ne permettent pas de 
dégager de claire corrélation avec les caractères moléculaires et morphométriques. Notre 
analyse morphométrique permet d'attribuer sans ambiguite le nom Rana Inmnocharis Gra- 
venhorst, 1829 à l'espèce à vaste répartition qui se trouve à Java, Sumatra et Bornéo Les deux 
taxa vivant en sympatric, nous donnons un nom nouveau, Fejemarya iskandari, à la popula- 
tion de Chianjur, Java. 
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Appendix 1 List of morphometric measurements of Fejervarya specimens used in this 


study. 
Abbreviation Measurement 
EL eye length 
EN distance from front of eye to nostril 
FFIF webbing measured from maximum incurvation between toes IV and V to tip 
of toe IV 
FLL length of arm from elbow to base of outer palmar tubercle 
FOL length of foot from proximal border of inner metatarsal tubercle to tip of toe 
IV 
FTL length of fourth toe from basal border of proximal subarticular tubercle to 
tip of toe IV 
HAL length of hand from base of outer palmar tubercle to up of finger II 
HL head length from mandibular articulation to tip of snout 
HW head width 
IBE distance between posterior borders of eyes _ 
IFE distance between anterior borders of eyes 
IMT length of inner metatarsal tubercle 
IN internarial distance 
ITL inner toe length from distal border of inner metatarsal tubercle to tip of toe 1 
MBE distance from mandibular articulation to posterior border of eye 
MFE distance from mandibular articulation to antenor border of eye 
MN distance from mandibular articulation to nostril x 
MTFF webbing measured from distal border of inner metatarsal tubercle to 
maximum incurvation between toes IV and V 
MTTF webbing measured from distal border of inner metatarsal tubercle to 
maximum incurvation between toes III and IV 
SVL |snout-vent length 
TFL length of finger from basal border of proximal subarticular tubercle to tip of 
finger HI 
TFTF webbing measured from maximum incurvation between toes ll and IV to 
tip of toe IV 
TL length of tibia from tibio-metatarsal articulation to knee 
maximum tympanum diameter 
tympanum-eye distance 


Source . MNHN, Paris 


Appendix 2 Allele frequencies, average expected heterozygosity (He) and polymorphism (P) at 15 presumptive enzyme loci an 9 samples of 
Feyervarya from Java (J), Kahmantan (K) and Sumatra (S) See text (Material and methods) for details on localmes Sample Ke belongs to 


the species Fejervarya cancrivora 
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Appendix 3 — Morphometrical data for adult specimens of Fejervarya samples from Sunda Islands Mean, 
standard deviation, minimum and maximum values are given for all samples 
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28 ALYTES 19(1) 
Appendix 3. — Continued 


Java 

Sidikalang Desa Lape Pontianak 

(Sumatra) Kalimantan) | (Kalumantan) | À ete 
1 female females 8 females 2 females Her 


1161473 108 = 104 1173397 116427 
FURY 111-121 107-109 109-120 14-118 
391404 893234 881 393 94 + 107 
AE 85-93 87-91 81-91 94 95 
2171577 2242379 2201779 151585 
ELL/SVL 212-225 222-227 207-230 221-229 
345 £1831 3482773 52452813 54771837 
PENY 519 563 543 -554 465- 558 534 - 500 
330 £1013 312 +096 3132 1846 336411 74 
Bei 315-338 331-133 273-337 128-144 
190 = 14 08 207499 201.511 87 208 + 4.94 
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39121267 316199 393 = 1807 385 + 187 
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